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Abstract 

We study theoretically the electron transport in a ID conductor adiabatically connected to a 
superconducting and normal metal leads. In the case of non-interacting electrons in the wire, we 
obtain full I-V curve. We show that ac voltage applied along with dc voltage modifies I-V curve 
resembling a regime of photon-assisted tunneling. We obtain in the limit of low voltages, that the 
electron-electron interaction affects electronic transport resulting in a decrease of conductance and 
features in the frequency dependence of the impedance. 
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I. INTRODUCTION 



It is well-known for a long time that a nonlinear conductivity arises in superconductor- 
normal metal contacts with multimode constriction, in case of both dirty constrictions 1 and 
constrictions without scattering^. This non-linearity is related to Andreev reflection. One 
can also expect nonlinear I-V curves in a junction of a superconductor with a quantum wire. 

In contrast to narrow ballistic quantum wires connected to wide normal metal reservoirs 
where the conductance is quantized in units of 2e 2 /h££, it was shown that in a wire with 
non-interacting electrons, which is connected to superconducting and normal leads Andreev 
reflection results in a double conductance quantum at zero temperature^. For the wire with 
normal leads electron-electron interaction does not modify I-V curves^-, although the simple 
Fermi liquid description with single-electron quasiparticles breaks down, and the Luttinger 
liquid (LL) picture that describes ID system of interacting electrons becomes an alternative 
to the Fermi liquid (for a review see Ref. b| and lldl). However, in a carbon atomic wire 
connected to a superconducting Al lead and to a normal Al lead density functional theory 
calculations^ predict reduction of the conductance due to electron-electron interaction. 

In this paper we study electronic transport in a quantum wire adiabatically connected to 
superconducting and normal leads. We derive boundary conditions for an adiabatic contact 
of a quantum wire with a superconducting lead. Then we calculate full I-V curves without 
taking into account electron-electron interaction in the quantum wire. We also consider static 
I-V curve in the case when the applied voltage contains dc and ac components simultaneously. 
The result resembles the regime of photon-assisted tunneling^ like in the systems with 
multimode constrictions. Finally, we study how the interaction affects conductance and 
obtain analytically contribution to the conductance due to electron-electron interaction. 

Below we set e, h and ^ to unity, restoring dimensional units in final expressions when 
necessary. 

II. BOUNDARY CONDITIONS 

We consider a long (in comparison with the Fermi length) ID conductor connected adi- 
abatically to superconducting (S) and normal (N) massive leads, see Fig. [TJ As the width 
of the quantum wire is of the order of the Fermi wavelength and small compared to super- 
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FIG. 1. System in study. S stands for superconducting lead, N denotes normal metal lead 

conducting correlation length, we can neglect suppression of the superconducting gap in the 
superconductor near the quantum wire due to the proximity effect. Using the quasi-classical 
approach one can justify the replacement of the adiabatically widening high- dimensional 
superconducting lead with an effective ID lead with a non-zero order parameter A 7^ at 
x < 0, as it was done in Ref.- 1 ^. Note that in case of normal metal leads the approach 
based on the replacement of the massive lead by an effective ID leacr^ and more detailed 
treatment of the boundary with a massive leao-^ yield similar results. 

At x > we describe interacting electrons by Tomonaga-Luttinger model (for a review 
gggSyxTj r-Q^ p era t or at x > can be represented as 

if.(x, t) = * St+ (x, t)e^ x + tt 8 ,_(x, t)e-^, (1) 

where s relates to spin, qp is the Fermi wave-number, and ^±(x) are field operators of 
right (left) moving electrons, and operators 4f±(x) slightly vary at distances of order q^ 1 . 

Boundary conditions for the contact of a quantum wire with a normal lead were derived 
in^. In the energy representation for an adiabatic contact these boundary conditions read 

* ai -{x = L,e)=R?{e), (2) 
Rf(e) = Ec s ,k2n6(e~e k ), (3) 

fc<0 

where c s> k are annihilation operators of the incident electrons with spin s and momentum k 
in the normal lead, e k is a spectrum of electrons in the normal lead. 

In the rest of this section we focus on the contact with a superconductor. In the super- 
conducting region field operators read 
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* t (x, £) = £ TfctM*) 6 "*** - 7^(x)e l£ ^, (4) 
*t(z, t) = £ 7jUfc(aOe fck * + 7fc t ^(^)e- i£fct , (5) 

where 7fc jS are annihilation operators of quasiparticles with momentum k and spin s, and 
u k (x), v k (x), £fc are corresponding solutions and an eigen energy of Bogolyubov equations 



(-I-I)"M +A "M = e "M> < 6 > 

A*m(»)- \-^-^) v ( x ) =m{x). (7) 



2m 2m r 

Solving Bogolyubov equations ©-(El) on a half-line x < 0, one can find a set of standard 
solutions 

u k (x) = (1/2 + ^/(2e))y 2 e^ 2 e ikx , 
v k (x) = (1/2 - ^ k /(2e)) 1 / 2 e~^ 2 e lkx J 

with an eigen-energy above the gap e k = a/|^| 2 + A 2 , where £ k = k 2 /(2m) — k 2 F /2m, x 
is a phase of the complex order parameter A. The solutions with k > kp and — kp < 
k < correspond to right-moving quasiparticles, while the other solutions correspond to 
left-moving quasiparticles. Beside of these solutions, there are two set of solutions with 
an eigen-energy below the gap, which decay into the depth of the superconductor. We 
parametrize these solutions by a decay parameter k 

_ J_ =F««/2+ix/2 ±ik F x+KX 

U>±,IK i/2 ' 

1 r ±a/2— i\/2 r ±ikpx+KX 



where = ^fivpK, e ia = {ivpn + e)/A. 

For the purpose of simplicity, below we take x — 0- 

Since the field operators 4f g (x) and their spatial derivatives continuously depend on x, 
we equate the operators in ([1]) and their spatial derivatives at x = with the corresponding 
operators in (S])-©. For the purpose of convenience we transform operators to the energy 



representation ^(e) = / ^>{t)e 1£t dt. Finally, after some algebra we relate fields 4f S}1 . at the 

— oo 

boundary of ID conductor with creation-annihilation operators of right-moving quasiparti- 
cles in the superconductor: 



Rf(-e) 



(8) 
(9) 



where £ R (e) = y/JjT+ iO) 2 — A 2 is analytic on the complex plane with a cut: (—A, A) 
and the right hand side is given by relations 



2^ 



R 



R 



.k>kp 



-e-A) ll,A £ + £ k) 



k=—kf 



Rf(-e) 



it 2£ R 
= 2tt- 



The correlation functions for these operators read: 



e-A) J%kS(e-e k ) + e(-e-A)J2^l,A £ + £ k) 

k=—kp k>kp 



R^ (s )R^(^E / 



RU-e')Rf\-e)) = 2 ^ l + ^ n F (e)27r5(e - s'), (10) 



where np(e) is the Fermi distribution. Note that if we take e = the boundary condi- 
tions (151)-® are reduced to those obtained in 13 . 



III. NON-INTERACTING CASE 

In order to obtain I-V curves for all voltages we assume first that electrons in quantum 
wire do not interact. We study the case when the voltage consisting of dc and ac components 
V{t) — Vq + V\ coswo^ is applied to the normal lead. The case of time-independent voltage 
can be obtained by taking V\ = 0. In the presence of time-dependent potential Heisenberg 

operators R£ (t) in the right hand of boundary conditions © should be replaced by R^(t) = 

t 

R^(t) exp(— i f V{r)dr). Thus, a correlation function in the time representation reads 
o 

iVo(tf — t) + i— - (sinwo^' — sincuo^) • 



(R^(t')R?(t)) = (R N \t')R N (t))e W 
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+00 , x 

After substituting fl3]) and using the expansion exp(i^smcoot) = Yl Jk (^) exp(— ikuot) 

k=—oo ^ ' 

we obtain the correlation function in the energy representation 



x J fc ( — J J m ( — ) 27r5(£: - e' - fcwo + mw ), (11) 



m,k 

where np{e) is the Fermi distribution. In the case of time-independent potential (Vi = 0) 
the correlation function reduces to 

G N (e, e') = —n F (e - V )2n5(e - e'). 
v F 

In the absence of interaction in the quantum wire operators iff sr (x) in the general ex- 
pression (JTJ read 

i> s>r (x, e) = *., r (0, e)e^, g(e) = — . (12) 

vf 

Using boundary conditions (j2J), (jEJ) - ® and expressing pairings of fermionic operators 
(W^) at x = in terms of correlators ( TTU1) . ( ITT1) and then transforming into the time 
representation we obtain the current j(t) = Vp Yl r (^tr(^)^s,r(t))- In the case of 

s=t4_,r=±l 

stationary potential the current reads 



Io(V) = G 



o 



6 (A — e) + 0(e - A)- 1 



tanh tanh as, 

2T 2T 



where Go — 2e 2 /h is the conductance quantum. In the limit of T = integration yields 

\V\-A), (13) 



2G 

I„(V) = 2G V6(A - \V\) + — - 



V 3 - (V 2 - A 2 ) 3 / 2 
{ —- 2 ] — + 2A sign (V) 



A 2 

The I-V curve is linear with double conductance quantum 2Gq at the voltages below the gap 
\V\ < A, and behaves sufficiently non-linearly at the voltages V <^ A and at high voltages 
V ^> A it is almost linear with a normal conductance quantum Go and an excess current 

lexc — jjGqA. 

In the limit T ^> A the I-V curve differs from the normal-state Ohm's law by the term 
that describes small correction to conductance at low voltages and saturates at high voltages 
giving the excess current 

Io(V) = G (V + ^Atanh^V 
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FIG. 2. Dependence of differential conductance on the applied voltage. Solid line corresponds 
to the dc voltage. Dashed line corresponds to the voltage consisting of dc component V and ac 
component Vicosuiot (the conductance was calculated for the case eV\ = Hujq = A) 

If an ac voltage is also applied then the averaged over time current is given by a Tien- 
Gordon-like^ formula 

+oo 

i{V)= h{v + khu )j 2 k 

fc=— oo 

The dependence of differential conductance on the applied voltage at T = is shown in 
Fig. El 

Originally Tien-Gordon formula was derived for tunnel contacts, however, as we see 
it is applicable in case of superconductor and normal leads connected by a conducting 
constriction as well. This occurs because the current in the constriction does not affect the 
electronic states in the lead, and the applied voltage is taken into account by the factors 
expz[ J dtV(t)} added to Green's functions in boundary conditions. As the current, in case of 
non-interacting electrons, linearly depends on the Green's functions in the leads, this yields 
the results similar to photon-assisted transport. In case of interacting electrons, when the 
Hamiltonian in the quantum wire is not bi-linear in electron operators, the Tien-Gordon-like 
relation is not expected. 

Up to now we considered the superconducting electrode as a natural superconductor, 
however, superconductivity can be induced in the two-dimensional electrode by the proximity 
effect. The proximity effect results in the superconducting pairing in the normal electrode 
described by the self-energy terms in the Green's function of the electrode^^. In this case 
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the energy gap A induced in the 2D lead by the superconducting metal is smaller than the 
gap in the superconductor Ag. In the limit of small tunneling rate T ~ * dk F ^ /\ s ^ w here 
t is tunneling matrix element, and d is the thickness of the electrode, A w T. It can be 
shown using the approach of Ref.— that the equations (IHl fTOi) can be used with substitutions 

e + 



e — > e 



irAs 



and £ R — > \/ e 2 — T 2 . Therefore, our results are applicable to the 



case of the proximity-induced superconductivity if we substitute A in previous equations for 
A « T. However, new features in the conductance are expected at voltage V ~ A5 ^> A 
where the non-linear part of the I-V curves (Tl3|) is small. The conductance at this voltages 
reads 



dl 
dV 



Gr 



i + n- 



r 2 (A| - v 2 



W 2 [T + y/A 2 s - (V + iO) 2 } 2 
When the voltage approaches to Ag the non-linear correction to the conductance (the second 
term) has a dip of the approximate width T 2 /As with a minimum value — GoT 2 /32V 2 at 
V w A s + T 2 /6A S . 



IV. INTERACTING CASE 

In this section we will study the effect of electron-electron interaction on the conductance 
when the applied voltage \V\ < A and V does not depend on time. We describe interacting 
electrons using Luttinger liquid (LL) model and bosonization technique 9 in terms of boson 
displacement field & u (x,t) and its conjugated momentum fl u (x,t) = d x Q u (x,t). Here v = 
p,a stands for charge and spin sectors correspondingly, and @„ is a field dual to & u . The 
standard Tomonaga-Luttinger(TL) Hamiltonian reads 

2 7>?_ / „ \ 2l 



l/=p,cr 



(14) 



Here the LL parameters K u playing the role of the stiffness coefficients of the elastic string 
described by Hamiltonian (}T4"|) , are related to the electron-electron interaction potential, and 
measure the strength of interaction between electrons. We consider spin-rotation invariant 
case in our study and short-range electron-electron interaction. Under these assumptions 
K a = 1, and K p = g ^~^ where g(q) is the Fourier-transformed interaction potential. 

Parameter K p = 1 corresponds to non-interacting case, and < K p < 1 for a repulsive 
interaction potential. 



In the non-interacting case considered in the previous section only the states below the 
superconducting gap A contribute to a current if the applied voltage \V\ < A. In the 
interacting case as distribution function does not display a sharp jump at the Fermi level, a 
contribution from the states above the gap becomes significant. We study the case of weak 
enough interaction 1 — K p <C 1 when this contribution can be considered as a perturbation. 

In the zero-approximation we neglect the contribution of states with energy \e\ > A. 
Therefore, the right-hand the of the boundary conditions (jHj)-© can be omitted, and since 
for \e\ < A the Bogolyubov amplitudes \u(e)\ = \v(e)\, the boundary conditions 
yield 

:ft + (i)t t)+ (i):=:t(i)*; r (t):, (15) 
: ^ j+ (t)*} i+ (t) :=: ^_(t)# tj _(t) : . (16) 

Here : A := A — (A)q denotes normal-ordering, (• • • )o stands for averaging over equilibrium 
state at V = 0. The boundary conditions (fI5]) - (fl6"]) can be bosonized using the relations 

:q?l r (t)y s>r := ^J- - s-^ + r * p + rs- 



2\f2-n 2V2ir 2y/2nv F 2^/2nv F 
resulting in the boundary conditions for displacement fields $j, and dual fields Q u 

d x $ p (x = 0) = 0, (17) 

dA{x = 0) = 0, (18) 

d t e p {x = o) = o, (19) 

d x e a (x = 0) = 0. (20) 
The boundary conditions for the normal contact were derived in^ (and earlier ire) 

- %d x Q p - dt$ p = P P , x = L, (21) 

- v F d x $ a - d& a = P a , x = L, (22) 
-v F d x Q p - d t Q p = -P p , x — L, (23) 
-v F d x e a - d t Q a = Pfj j x = L, (24) 

where P p ^(t) = 2| [: Ftf r \t)Ftf r (t)±R?\t)R?(t) :] , (P p ) = -V2V and operators B% are 
defined in (13]). 



The TL Hamiltonian ( TH|) yields equations of motion for the fields Q u 

- $>d 2 A = 0, (25) 

V 

- % d l®« = °- (26) 

V 

It is convenient to decompose the displacement field $ p and the dual field P in the sum 
of thermodynamically averaged c-number components <p p = ($ p ), p = (Q p ) and fluctuating 
components <p p — <fr p — <p p , § p = Q p — 9 P . Solving wave equation fl25|) subject to boundary 
conditions 0Tj) - ([TSl) . (l2"Tj) - (l2"2l . we obtain a zero-approximation solution for <f> p : 

<\> p = y/2Vt, (27) 

^(x = 0,u)= , . ,„ r/^r^ K n ^ ( 28 ) 

' ' uj(sm{K p ujL/vp)+iKpCos KpU)L/vp) \ ' 

The current can be calculated as j = \/2dt(f)p/iT = 2GqV. Thus, zero-approximation solution 
for interacting system at low and constant voltages yields the same conductance as in non- 
interacting case. However, if ac voltage with a frequency u is applied V = V\ coscjo^, then 
{P p (t)) = — V2Vi cosu t, and the averaged solution for the wave equation (1251) subject to 
boundary conditions flT A - fllS]! . (l2T]) - (j2"2]) reads 



<pp{x,t) = . r7 — r COSK P U X V F + C.C. 

Uq {sm{KpU) Q L/vF) + iK p cos K p qu L/vf) 
Thus, even in the zero approximation which neglects contributions from energies |e| > A, 

the electron-electron interaction affects an impedance Z(uS) of the system 



Z(u) 



cos K p uqL/vf — iK p 1 sin K p uqL/vf 



2Gq cos 2 uL/2vf 

It is worth noting that if an ac voltage is applied to the right electrode, since there is an 
average potential Uq = V^/2 the wire is being charged and discharged periodically, and the 
current flowing into and out of the wire is affected by the interaction 
/„ = /„(0) - L(L) = 2G„U Q . ^.L/2 VF 

Correlation functions and commutators for fluctuating parts <p p , d a in case of T = and 
a long wire (L ^> A;^ 1 ) read 

({^(0),^(r)}) = 2(^(0),^(0))-2ir p lnA|r|, (30) 

[<&>(°)> Vpi. T )} = 2iK p sign r, (31) 
'^(0),^(r)}) = 2(^(0), ^(0))-21nA|r|, (32) 
[<p a (0), <p ff (r)] = 2i sign r, (33) 
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where A is cut-off parameter and is of order of £p, and |r| ^> A -1 . Note, that since the TL 
Hamiltonian (j!4p is quadratic in $ v and the boundary conditions ffT71) - fT22l) are linear, the 
fluctuations of the displacement fields $ v are Gaussian in zero approximation. Fermionic 
fields ty sr are related to the field Q v by equation 



* Sir = [W^-exp 



(34) 



where i7 rs are unitary Klein factors which decrease the number of fermions by one and 
commute with bosonic fields $ y and 6^. Using f[3"4"|) . relations f[3"Uj) - f[3"3l and the identity 
e A e B = e A + B e l A < B }/ 2 we obtain fermionic correlation functions at the boundary x = in time 
and energy representations 

V(#U + r)) = (n, + (m, + (t + r)) = ^e^ ^gy 1 , (35) 
= (*| f+ (0*t. + (e)) = - ^) A(1 X P)/2 ^fe^ ^W. ( 36 ) 

+oo A/|e| 

F(^) = sinvr^- / ^^^-cosvr^- / ^^dx. (37) 



In the rest of this section we focus on the case of dc voltage and take into account the 
contribution to the current due to the states above the superconducting gap considering 
this contribution as a perturbation. The fermionic field ^ at the boundary x = can be 
represented as sum of the solution in zero approximation and a perturbation: ^ = x|/( )+^( 1 ). 
We multiply fermionic boundary conditions (JHJ) — (EJ) by the complex conjugate, and take into 
account that /^f^pR^ = since operators Rf and R^j are independent and the zero- 
approximation solution can be expressed in terms of Rf using ( l28l) . Finally, we obtain 
boundary conditions for the first-order terms 

^ (o)\ (38) 
: *\ + (e')* t+ (e) :) J + (: R S ^(e')R^(e) :), 

: *\_(e')%„(e) :) W - (^(^U^ = - A)| ((: % + (e')*\ + (e) ^ + 
: *\_(e')%4s) :} ° J - <: J2f (e')i?f ( e ) :). 

If the potential is applied to the normal lead then R^ (e')Rg (e) :y = 0. Using the 
values of fermionic correlators in zero approximation f )36|) we transform fl38|) - fl39l) to time 
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representation at low voltages V < A 



( 




where A(K p ) = J 




dx. Thus, the contribution from the states above the gap can 



be considered as small only if interaction is weak enough (1 — K p ) ln(A/A) -C 1. In this case 
F(K p ) w 1/4 and A(K p ) « ln(A/A). The boundary conditions for the first-order terms f )40|) 
in the bosonized form for K p w 1 read 



Thus, we find an averag ed solution <f>f> = ^2V{1 - K p ) ln(A/A) • {x/v F - t). Then the 
contribution to the conductance in the first order of perturbation theory reads 



V. CONCLUSIONS 

We have studied electronic transport in ID conductor attached to a superconducting 
lead and to a normal metal lead. If electron-electron interaction in the ID conductor is 
neglected, I-V curve is linear at low voltages V < A/e and is sufficiently non-linear at 
V > A/e. In the limit of high voltages V 3> A/e I-V curve is almost linear with an excess 
current |G A/e If the applied voltage contains ac component then the static I-V curve is 
modified resembling the regime of photon assisted tunneling^. We show analytically that 
in the case of a weak interaction in the quantum wire, the interaction results in features in 
the frequency dependence of the impedance and reduces the dc conductance. 




(41) 
(42) 




(43) 
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